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Abstract. The Atacama Cosmology Telescope has measured the angular 
microwave fluctuations to arcminute scales at frequencies of 148 and 218 
seasons of data. At small scales the fluctuations in the primordial Cosmic 
ground (CMB) become increasingly obscured by extragalactic foregounds and secondary 
CMB signals. We present results from a nine-parameter model describing these secondary 
effects, including the thermal and kinematic Sunyaev-Zel’dovich (tSZ and kSZ) power; the 
clustered and Poisson-like power from Cosmic Infrared Background (CIB) sources, and their 
frequency scaling; the tSZ-CIB correlation coefficient; the extragalactic radio source power; 
and thermal dust emission from Galactic cirrus in two different regions of the sky. In order to 
extract cosmological parameters, we describe a likelihood function for the ACT data, fitting 
this model to the multi-frequency spectra in the multipole range 500 < l < 10000. We extend 
the likelihood to include spectra from the South Pole Telescope at frequencies of 95, 150, and 
220 GHz. Accounting for different radio source levels and Galactic cirrus emission, the same 
model provides an excellent fit to both datasets simultaneously, with x 2 /dof= 675/697 for 
ACT, and 96/107 for SPT. We then use the multi- frequency likelihood to estimate the CMB 
power spectrum from ACT in bandpowers, marginalizing over the secondary parameters. 
This provides a simplified ‘CMB-only’ likelihood in the range 500 < i < 3500 for use in 
cosmological parameter estimation. 

Keywords: cosmological parameters from CMBR, CMBR experiments, Sunyaev-Zeldovich 
effect 


power spectra of 
GHz, from three 
Microwave Back- 


ArXiv ePrint: 1301.0776 


J CAP 0 7 (2013) 025 


Contents 


1 Introduction 1 

2 Model for the mm-wave sky 2 

2.1 Thermal Sunyaev-Zel’dovich 5 

2.2 Kinematic Sunyaev-Zel’dovich 5 

2.3 Cosmic infrared background 6 

2.4 tSZ-CIB cross-correlation 7 

2.5 Radio point sources 8 

2.6 Residual Galactic cirrus 9 

3 Full likelihood from small-scale data 9 

3.1 Likelihood from the ACT data 9 

3.1.1 Calibration and beam uncertainty 11 

3.1.2 Secondary model parameters 11 

3.2 Combining with SPT data 11 

3.3 Multi-frequency likelihood prescription 13 

3.3.1 Combining with WMAP 13 

4 Tests of the multi-frequency likelihood 13 

4.1 ACT data 15 

4.2 Combination with SPT 16 

4.3 Tests of the likelihood 17 

5 CMB-only likelihood 19 

5.1 Method: bandpowers via Gibbs sampling 20 

5.1.1 Combining spectra from different regions 22 

5.1.2 Calibration factors 22 

5.2 Marginalized CMB bandpowers 22 

5.3 The CMB-only likelihood 23 

6 Summary 24 

A Tests of the CMB likelihood 25 


1 Introduction 

Measurements of the Cosmic Microwave Background (CMB) have played a central role in 
constraining cosmological models. Anisotropies measured over the whole sky by WMAP 
have provided evidence for a flat universe described by just six cosmological parameters. 
The measurement of the Sachs- Wolfe plateau in the power spectrum, and three acoustic 
peaks, have led to constraints on ACDM parameters to percent-level accuracy ([38, 40]). 
The Silk damping tail of the power spectrum provides a wealth of additional information 
about the physics of the early universe, encoded in its shape, and in the positions and 
heights of the higher-order acoustic peaks ([64]). Extracting information from these angular 
scales is complicated by the presence of additional power from extragalactic point sources, 
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Dataset 

Frequency 1 

GHz 

Reference 

Area 

sq degrees 

•^min 

^max 

Sc 2 
mJy 

t'tsz 3 

GHz 

I4tad 

GHz 

i'CIB 

GHz 

ACT 

148 

[14] 

590 4 

500 

10000 

15 

146.9 

147.6 

149.7 


218 



1500 

10000 


220.2 

217.6 

219.6 

SPT-low 

150 

[36] 

790 

650 

2000 

50 

152.9 

150.2 

153.8 

SPT-higli 

95 

[54] 

800 

2000 

9400 


97.6 

95.3 

97.9 


150 



2000 

9400 

6.4 

152.9 

150.2 

153.8 


220 



2000 

9400 


218.1 

214.1 

219.6 


Table 1 . Small-scale CMB datasets 


emission from the Galaxy, and secondary anisotropies due to the thermal and kinematic 
Sunyaev-Zel’dovich effects ([67]). 

The Atacama Cosmology Telescope (ACT) mapped the mm- wave sky with arcminute 
resolution from 2007 to 2010 in two distinct areas. About 600 square degrees were used 
to compute the angular power spectrum. Power spectra and cosmological results using the 
1-year data, from the 2008 observing season, were presented in [13, 22] and [16]. During 
roughly the same period, the South Pole Telescope also mapped the microwave sky, and 
presented cosmological results in [36, 43, 54, 60] and [66]. 

In this paper we describe a method to fit multi-frequency power spectra from the ACT 
data simultaneously for CMB, foreground, and SZ parameters, following a similar approach 
to analyses in [16] and [54]. We describe the likelihood constructed for the 3-year ACT 
dataset, using data from the 2008-2010 observing seasons, and show how it can be used in 
combination with data from SPT in a self-consistent way. Using this likelihood from ACT, 
we then construct a simpler CMB-only likelihood, estimating CMB bandpowers marginalized 
over the SZ and foreground parameters. 

This is one of a set of papers on the 3-year ACT data; [14] present the angular power 
spectra, and [61] use the likelihoods presented here to estimate cosmological parameters. We 
begin in section 2 by describing the model for the mm-wave emission. In section 3 the full 
likelihood of the ACT data is described, including the combination with data from SPT and 
WMAP 7-year data. In section 4 we show the small-scale model fit to the multi-frequency 
data. In section 5 we describe the compressed CMB-only likelihood, concluding in section 6. 

2 Model for the mm-wave sky 

Sky maps of mm-wave fluctuations at arcminute resolution include components emitting at 
low redshift, in addition to the primordial CMB signal and secondary CMB effects, e.g., [16, 
22, 36, 63]. The power spectra from the complete ACT dataset, reported in [14], are shown 
in figure 1, focusing on angular scales of interest for the primordial CMB signal. At scales 
smaller than a few arcminutes (t ~ 1500) the secondary signal, which we define as the sum of 

X A11 cross-spectra between channels are used in the likelihood. 

2 F1ux cut imposed on map by its point source mask. 

3 Effective band-centers from ACT are from [68], given for tSZ, radio sources, and CIB sources. 

4 This area includes the ACT-E region at dec = 0° (300 deg 2 ), and the ACT-S region at dec = —55° 
(290 deg 2 ). 
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2000 



£ 

Figure 1 . Summary of small-scale mm- wave data measured by the Atacama Cosmology Tele- 
scope ([14]) and the South Pole Telescope ([36, 54]), in the angular range used for measuring the 
damping tail of the CMB. The ACT and SPT data are independently calibrated to WMAP. The 
vertical axis is £ A Ce instead of the conventional £ 2 Ct to highlight the features at these angular scales. 
The primary CMB signal corresponding to the best-fitting ACDM model ([61]) is indicated (dashed), 
together with the total signal at 148 GHz (red, lower solid curve) and 217 GHz (black, upper solid 
curve), including secondary effects from SZ and foregrounds. Modeling the secondary contributions 
from SZ and foregrounds is vital to allow extraction of the primordial signal at small scales. 



Figure 2. Regions of the sky used for ACT power spectra (dashed, [14]) in the Equatorial plane 
(ACT-E, 300 deg 2 ), and at -55° declination (ACT-S, 292 deg 2 ). The 800 deg 2 used for SPT power 
spectra (dotted, [36, 54]) is indicated, with 54 deg 2 overlap with ACT-S. The color scales with Galactic 
cirrus intensity ([21]). 


foregrounds and SZ effects, becomes significant compared to the CMB. We want to extract 
the primary CMB signal, but, since there are more foreground components than frequency 
channels, information about both the frequency and scale dependence of the foregrounds is 
required to separate the signals. In this section we describe a model to fit the power spectrum 
of these fluctuations over the frequency range 90 ^ v ^ 250 GHz probed by ACT, SPT, and 
other CMB experiments including the Planck satellite ([50]). We follow a similar approach 
to [16, 54, 63]. 
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For frequency v and direction n we model the signal in the maps as 


A T{y, n) = AT CMB (n) + A T se >, n), (2.1) 

where AT CMB (n) are the lensed CMB fluctuations, which are independent of frequency in 
thermodynamic units. The secondary signal, AT sec (y n), is dominated by the sum of tSZ 
and kSZ components, emission from dusty infrared galaxies and radio galaxies, and dust 
emission from Galactic cirrus, all of which are functions of frequency. 

The cross-correlation power spectra between frequency i \ and v» are calculated as 

cj ® (r^i/OTKi/j)) , (2.2) 


where Tf is the Fourier transform of T(n) in the flat-sky approximation. The theoretical 
cross-spectrum S)* 1-1 ' = £(£ + l)C^ h,1 Y27T is modeled as 

£f’ ij ^ MB +sr ,ij > ( 2 - 3 ) 

where Bf MB is the lensed primary CMB power spectrum. In this analysis we model the 
secondary spectra as 

0f c,ij = Sf Z;ij + Bf : Z,ij + £° IB_P ’ ij + £° IB_C ’ ij 

+fif Z-CIB,ij + B rad,ij + ^Gal.ij ^ ( 2 . 4 ) 


with contributions from the tSZ and kSZ effects; dusty galaxies that form part of the Cosmic 
Infrared Background, both Poisson-like (CIB-P) and clustered (CIB-C); the cross-correlation 
between the tSZ effect and the CIB (tSZ-CIB); radio galaxies (rad); and dust emission from 
Galactic cirrus (Gal). We assume that all other cross-spectra can be neglected. Measure- 
ments by WMAP and other CMB experiments show that other Galactic emission, including 
synchrotron and free-free emission, is negligible in the v £ 90 GHz frequency range at these 
small scales and locations, e.g., [25]. The cross-correlation of radio sources and both tSZ 
and CIB sources is also expected to be small, e.g., [55], who find a correlation of only a few 
per cent in simulations. Since the kSZ signal consists of positive and negative fluctuations, 
depending on the line-of-sight motion of the electrons that source the signal, the two-point 
correlation function with other signals should average to zero. 

The majority of these secondary spectra will be common to all regions of the sky and 
to different experiments. The power in the residual radio point sources is expected, however, 
to vary among data sets due to the removal of bright sources. For example, for a Poisson 
distribution of sources with differential number counts scaling as dN/dS oc S~ 2 , the Poisson 
power will be 

C, = l S>- dS , (2.5, 

i.e., Cf oc Sjnax, where Sjnax is the flux of the brightest sources in the map at a given frequency. 
Deeper surveys, with lower noise per pixel, are able to detect and mask out dimmer sources, 
so S’ma.x will be lower, leading to a lower residual power. Radio sources have a shallower 
dN/dS slope than CIB galaxies, so imposing a flux cut of, e.g., Sjnax = 15 mJy removes a 
significant amount of radio power, but little CIB power. 

The Galactic foreground power is also expected to vary between regions on the sky. 
The two regions mapped by ACT (Equatorial, ACT-E, and South, ACT-S) are shown in 
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figure 2 and summarized in table 1, together with data from SPT. The temperature scale 
of the Galactic cirrus map from [21] is shown for comparison. A higher level of emission is 
expected in some regions of the ACT-E region. 

In the rest of this section we describe how each of the components in eq. (2.4) are 
modeled. To allow for comparisons between experiments, we normalize the power spectra 
at a pivot frequency of t'o = 150 GHz and scale — 3000. This differs slightly from the 
convention used in previous analyses of the ACT data ([16, 22]). In each case we describe the 
parameterization used in the fiducial model; in a later section we consider possible extensions 
or modifications. 


2.1 Thermal Sunyaev-Zel’dovich 


Our model for the power from thermal SZ fluctuations is given by 


r? tSZ,ij 


= atsz 


/(^)/(^) 

/> o) 


B 


tsz 
0 / > 


( 2 . 6 ) 


where Bq SZ is a template power spectrum corresponding to the predicted tSZ emission at uq 
for a model with amplitude of matter fluctuations erg — 0.8, normalized to 1 p, K 2 at £o — 3000, 
and atsz is a free parameter describing its amplitude. An example is shown in figure 3. The 
factor f(u) — xcoth(x/2) — 4, for x — hu /ksTcMB, scales the expected tSZ emission to 
thermodynamic units at v, the effective band-center for the tSZ, given in table 1. We ignore 
relativistic corrections, e.g., [34], since the low-mass clusters that dominate the spectra are 
well approximated by the non-relativistic formula. This atsz normalization differs from that 
used in [16]. The present choice has the advantage of reducing the dependence on the choice 
of template, since the main difference between various templates is their amplitude. This 
means one expects to find the same constraint on atsz regardless of template, and a constraint 
on the SZ power can be converted back into a model-dependent constraint on erg. 

The template we adopt is derived from recent hydrodynamic simulations described in [4] . 
The simulations include the effects of radiative cooling, star formation, and feedback from 
AGN and supernovae. The predictions are consistent with SZ measurements from both SPT 
and ACT, e.g., [43], and the shape is shown in figure 3. For the model with as — 0.8, the 
predicted spectrum reported in [4] has amplitude atsz — 5.6 ± 0.9, with standard deviation 
estimated from ten simulations. 

Numerous other authors have also predicted the tSZ spectrum from independent sim- 
ulations and analytical models, e.g., [4, 19, 37, 55, 57, 59, 69], and the expected amplitude 
for fixed cosmological model varies depending on the astrophysical modeling of the clusters. 
However, the template shape is broadly consistent among models, and the data are not yet 
sensitive to shape difference, so we do not include a shape uncertainty. We do not mask 
clusters, and expect the total SZ power to be the same for ACT and SPT. 


2.2 Kinematic Sunyaev-Zel’dovich 

The kSZ power is expected to have contributions arising from fluctuations in the electron 
density ([49]), and in the ionization fraction, e.g., [27, 33, 45], as well as from the motion of 
galaxy clusters at later times. We model the power as 

Bf ZM = a kS zB 0 kS ; Z , (2.7) 

where B^ sz is a template spectrum for the predicted blackbody kSZ emission for a model with 
erg = 0.8, normalized to 1 /iK 2 at £q — 3000. The parameter a k sz describes its normalization. 
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Figure 3. Template power spectra for the thermal and kinetic Sunyaev-Zel’dovich effects (tSZ and 
kSZ, [4, 5]), clustered CIB sources scaling as £ 08 (CIB, [2], the cross-correlation between tSZ and 
CIB (tSZ-CIB, negative at 150 GHz, [3], and Galactic cirrus ([47]). They are normalized at £ = 3000 
and 150 GHz, and the tSZ-CIB is shown for a perfectly correlated signal. Poisson CIB and radio 
source power (not shown) scale as i 2 . 

We use a template that assumes a model with instantaneous reionization, described in [5]. 
This is derived from the same hydrodynamic simulations as the tSZ spectra in section 2.1, 
and is shown in figure 3. The predicted amplitude from the simulations is aksz = 1-5 for 
homogeneous reionization at z — 10 in a <78 — 0.8 cosmology. This is a quarter of the 
expected tSZ power. The corresponding kSZ template for the ‘nonthermal20’ model in [69] 
has a similar amplitude and shape, as does the [58] ‘CSF’ model, and the [7] model. The 
power is expected to scale as roughly Cg" 5-5 ([58, 69]). 

Reionization of the universe is not expected to be instantaneous, as was assumed by [5]. 
The shape and amplitude of the kSZ power from patchy reionization is far less certain, with 
simulations predicting a signal at least as large as the homogeneous signal, e.g., [6, 46, 73]. 
This gives a total expected signal of ~ 3 to 5 /i K 2 for simple reionization models at i — 3000, 
comparable to the tSZ at 150 GHz. The dominant effect of patchy reionization on the power 
spectrum at scales probed by ACT is to alter the amplitude, depending on both the midpoint 
and duration of reionization, e.g., [6]. We test a modified shape in section 4.3, but do not 
include additional shape uncertainty in the template in the basic model. 


2.3 Cosmic infrared background 

Thermal dust emission from high redshift star-forming galaxies, part of the Cosmic Infrared 
Background (CIB), is emitted in the rest- frame far infrared and redshifted into the mm- wave 
range, e.g., [31, 53]. Clustering of these galaxies has been detected statistically in mm- wave 
maps at CMB frequency ([16, 29, 30, 52, 54, 60]), as well as in the sub-mm, e.g., [39, 71]. 
Following the analyses in [2, 16, 54], the power from these galaxies is modeled as the sum of 
a Poisson and clustered component, given by 


tjCIB— P ,ij 



'Vvi,Pp)Vvj,/3 P y 

Vo) 

M 2 (z ' o,/3 P ) 


^K 2 


( 2 . 8 ) 
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for the Poisson part, and 


t^CIB— C,ij 


( £\ 2 ~ n 

' H{vi,Pc)li{vj,Pc)' 

Vo) 

0,Pc) 


^K 2 


(2.9) 


for the clustered part. Here, n is a power law index, and the frequency scaling of each 
component is given by a modified blackbody, 


v(v,P) = B v {T d )g{y ), 


( 2 . 10 ) 


with emissivity indices f} p and /3 C for the Poisson and clustered dust terms respectively. The 
function B v {Td) is the Planck function at frequency v for effective dust temperature T<j, and 
the function g(v) — (dB v (T) / dT)~ l |t cmb converts from flux to thermodynamic units. The 
parameters a v and a c normalize the two components at £o and v o, and different frequencies 
channels are assumed to be perfectly correlated. 

The frequency dependence we adopt is an approximation to a sum of modified black- 
bodies at different redshifts, so this emissivity and temperature are only effective properties 
of the dust. Following [2] we fix the dust effective temperature to — 9.7 K. We also assume 
Pp — Pc in the basic model. 

The power-law angular scaling of the clustered term, with increasing £ 2 Ci power at 
small scales, is shown in figure 3, and approximates the shape of the non-linear power spec- 
trum, which includes contributions from pairs of galaxies in the same dark matter halo, and 
between galaxies in different halos. [2] find that a power law in £ provides a good fit to 
small-scale power spectra from Planck and the Balloon-borne Large- Aperture Submillimeter 
Telescope (BLAST), and from cross-correlating ACT and BLAST maps. This is consistent 
with observations of the correlation function from high-redshift Lyman break galaxies, as 
well as local galaxies, e.g., [12, 23]. We fix the power-law index to n — 1.2 in the fiducial case 
(£ 0 - 8 ), in close agreement with the estimate of n — 1.25 ± 0.06 in [2], Both are in agreement 
with galaxy correlation functions. 

We do not expect the CIB power to vary significantly between the ACT and SPT maps, 
despite the different flux cuts applied to remove sources. Using the model in [1], the predicted 
effect of source masking on the CIB power is only at the per cent level. 


2.4 tSZ-CIB cross-correlation 

Some spatial correlation is expected between clusters that contribute to the tSZ, and CIB 
galaxies, since both trace the matter density field. The higher redshift and lower mass groups 
that make an important contribution to the tSZ signal ([6, 37, 69]) are also likely to host 
dusty galaxies. [3] model this correlation, and predict the scale and frequency dependence 
of its angular power spectrum. For mm-wave spectra at £ > 2000, a correlation of ~ 10 to 
30% in power is predicted, with uncertainty dominated by uncertainties in the halo mass and 
redshift distribution of the CIB. A significant fraction of the CIB power on small scales is 
due to pairs of galaxies occupying group and cluster-mass halos (the ‘one-halo’ term of the 
halo model). These same halos are responsible for the tSZ power. A tSZ-CIB correlation 
in units of power of tens of per cent is therefore possible on small angular scales even if the 
overall fraction of CIB emission associated with massive halos is small. 

The tSZ-CIB power is negative at 150 GHz, and can partially cancel power from the 
kSZ effect as it does not vary significantly with frequency over the range probed by ACT and 
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SPT. As a result, neglecting this component can lead to artificially tight constraints on the 
kSZ power ([46, 73]). Following [2] we model the spectrum as 



( 2 . 11 ) 


J.OI7 piR • . o 

where B 0 ^ is the predicted correlation spectrum shape, normalized to 1 gK z at and 
shown in figure 3. The free parameter £ is the correlation coefficient. The Poisson CIB 


.tSZ-CIB 


parameter a p is not included in eq. (2.11), unlike in [60, 73], as the sources that dominate the 
CIB Poisson power in the mm-wave bands are unlikely to have significant redshift-overlap 
with the tSZ clusters ([3]). Assuming the same modified blackbody scaling for the CIB as 
in section 2.3, and tSZ frequency scaling f{v), the frequency scaling of the cross-spectra in 
thermodynamic units is then 


for pivot scale vq. 

Since the correlation coefficient is poorly constrained by ACT, we impose a uniform 
prior of 0 < £ < 0.2 in the basic model; the effect of widening the range to e.g., £ < 0.5, 
corresponding to the maximum allowed correlation in the models explored by [3] , is discussed 
in [61]. Due to the correlation between £ and the kSZ power, broadening the limit on £ 
increases the upper limit on the kSZ power, but does not affect cosmological results. 

2.5 Radio point sources 

The radio sources at ACT frequencies are not expected to be significantly clustered, see 
e.g., [30, 56], and to good approximation their power should be perfectly correlated between 
neighbouring frequencies, consistent with simulations in [55] and valid for sources with the 
same spectral indices. As in [16] and [54], we model the residual power after masking bright 
sources as Poisson scale-free power, with 


in thermodyamic units, where g{v) converts from flux units as for the CIB sources. The 
amplitude a s is normalized at vq and £q. Measurements of bright sources from ACT and 
SPT give an estimate for the spectral index in flux units of typically a s — —0.5 ([44, 70]). 
Assuming it holds at fainter fluxes, we fix a s — —0.5 in the fiducial model. 

Bright source counts can also be used to predict a s by extrapolating to fainter fluxes 
using a model for the number of sources as a function of flux. This was done in [44] for ACT, 
and for SPT in [36, 54]. Using point sources measured from the full ACT dataset, [26] now 
predict a residual power a s — 2.9 ± 0.4 after masking sources brighter than 15 mJy in both 
ACT regions (the level used to construct the mask for our maps in this analysis), where the 
catalog is estimated to be complete. We impose this as a Gaussian prior on the power at 
150 GHz. For comparison, the estimated power in the SPT power spectra after masking to a 
flux level of 6 mJy for SPT-high is a s — 1.3 ± 0.2 ([54]), which we also impose as a Gaussian 
prior. 


f'iyij ) = f(yi)n{vj,Pc) + 


( 2 . 12 ) 



(2.13) 
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2.6 Residual Galactic cirrus 


The Galactic emission is spatially varying, and [14] show that dust emission can contribute 
significantly to the power spectra, particularly in our Equatorial region. As reported in [14], 
we apply a mask to regions of high dust emission before computing the power spectrum, 
using measurements at 100 /. im from IRIS ([47]). 

We then marginalize over a residual Galactic cirrus component using a power-law 
template 


.-iGal.ij 



ViVj\ 

) 


09 


' g{vi)g{vjY 
. g 2 (vo) . 


mK 2 , 


(2.14) 


with amplitude a g , frequency index /3 g — 3.8, and angular scaling n g = —0.7. This angu- 
lar scaling is estimated from the 100 fim IRIS dust maps ([47]). The frequency scaling is 
estimated by correlating the IRIS dust maps with ACT ([14]), and is consistent with early 
results from Planck ([51]). Using the correlation coefficients estimated in [14], we impose 
priors of a ge — 0.8 ± 0.2, and a gs — 0.4 ± 0.2 in the ACT-E and ACT-S spectra respectively. 

For the SPT data, a small Galactic cirrus residual is also expected. In our basic model 
we follow the treatment in [54], fixing the Galactic cirrus power to £>3000 — 0.16, 0.21, and 
2.19 fxK 2 in thermodynamic units at 95, 150, and 220 GHz, with scale dependence oc £~ 1/2 . 
However, this model has a steeper angular power law than in our ACT model, and a shallower 
frequency scaling (/3 m 3.6 between 150-220 GHz). For consistency we therefore test the effect 
of adopting the ACT model instead, using £~ 0 - 5 * 7 and f3 g = 3.8, and a prior of a g — 0.4 ± 0.2. 
We find no effect on parameters and no change in the goodness of fit. 


3 Full likelihood from small-scale data 


In this section we describe the multi-frequency likelihood used to model the ACT data, and 
show how we extend it to include other small-scale datasets, in particular data from SPT. 


3.1 Likelihood from the ACT data 

The data from [14] describe the two ACT regions separately (ACT-E and ACT-S); and consist 
of multi-season and multi- frequency spectra, with an associated covariance matrix. ’ They 
are derived from ACT maps obtained using the method described in [18]. The likelihood 
is Gaussian-distributed to good approximation. To construct the likelihood for each region, 
given some model spectra we compute bandpower theoretical spectra using C^’ 1 ^ — 

where w l A is the bandpower window function in band b for cross-spectrum ij , 
described in [14]. 

The likelihood, «£?, of the data for each ACT region separately is given by 

- 21n2z? » ( Cf - C fe ) T S- 1 (C'6 h - C b ) + IndetS, (3.1) 


where E is the bandpower covariance matrix. Each of the model and data vectors C [ ) h and 
Cb contain three sets of spectra, 


_ rr ,148,148 .-,148,218 ,-,218,2181 
! ) L h l 


c b = \c;; 


5 In this paper we use the original version ‘vl’ of the spectra, unless stated. Since original release, the 

binning and the beams have been slightly refined, generating ‘v2’ of the spectra. This is described in [14], 

and the ‘v2’ spectra are released publicly; cosmological effects are negligible as described in [61]. 
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Figure 4. (Top) Power spectra measured by ACT ([14]) at 148 and 218 GHz, and their cross- 
spectrum, coadded over ACT-E and ACT-S. We show the primary (lensed CMB in dotted black line) 
and secondary contributions (dotted lines) to the best-fitting model. (Bottom) Residual power in the 
ACT cross-frequency spectra, after subtracting the best-fitting model, at 148 (left), 148x218 (center), 
and 218 GHz (right). The errors at small scales are correlated due to beam uncertainty. The model 
is a good fit simultaneously to ACT-E and ACT-S, with no sigificant residual features. 


for ACT-E and ACT-S separately, and each spectra set CV itself contains spectra for each 
cross-season. There are two seasons used for ACT-E (3 cross-season spectra), and three for 
ACT-S (6 cross-season spectra). The total likelihood is given by 

— 21nj£?ACT = — 2hl22?ACT-E — 2hi2z?ACT-S- 
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3.1.1 Calibration and beam uncertainty 

The data power spectra are calibrated, but have uncertainties. We therefore include a cali- 
bration parameter yi, for each map i, that scales the estimated data power spectra as 

Cf > y iyj cl 

and the elements of the bandpower covariance matrix as 

^bb' iyiVof^bb" 

To account for both ACT regions, we include four calibration parameters: y \ e , yi e for ACT-E 
at 148 GHz and 218 GHz, and y\ s , y 2 s for ACT-S. 

[14] calibrate the 148 GHz maps using WMAP, following the method in [28], at an 
effective £ — 700, resulting in a 2% map calibration error in CMB temperature units. We 
impose this as a Gaussian prior, with y\ e , y\ s — 1.00±0.02. The 218 GHz maps are calibrated 
relative to 148 GHz, at an effective £ = 1500. The 218 GHz calibration is constrained by the 
cross-spectrum, so no prior is imposed on y^ e and y^s- Within each frequency, the individual 
seasons are calibrated to each other; the inter-season calibration error is absorbed into the 
single overall calibration uncertainty. 

Uncertainties in the measured beam window functions for ACT at 148 GHz are between 
0.7 and 0.4%, and at 218GHz between 1.5 and 0.7%. We incorporate uncertainties in the 
measured beams by including them directly in the covariance matrix for the spectra, described 
in [14]. This technique assumes a fiducial model for the power spectra but is insensitive to 
its exact form. 

3.1.2 Secondary model parameters 

Our model described in section 2 has nine free secondary parameters for ACT in the basic 
case: atsz and atsz describing the SZ emission, a p , a c and /3 C describing the CIB power, 
a s describing the radio power, £ describing the tSZ-CIB cross-correlation, and a ge and a gs 
describing the Galactic cirrus emission. The latter four have strong priors imposed, as de- 
scribed in section 2: a s — 2.9 ± 0.4, 0 < £ < 0.2, a ge — 0.8 ± 0.2, and a gs — 0.4 ± 0.2. In 
addition to the nine model parameters, there are four calibration parameters for ACT. In 
section 4.3 we investigate how additional, or fewer, parameters affect the fit of the model to 
the data. To compute the model requires an effective frequency for each component; we use 
the band-centers for SZ, radio, and dusty sources given in table 1 ([68]). 

3.2 Combining with SPT data 

The South Pole Telescope observed the sky from 2007-10. Spectra are reported in [36] for 
angular scales 650 < £ < 3000 at 150 GHz, and in [54] for angular scales 2000 < £ < 9400 
at 95, 150 and 220 GHz. These observations are summarized in table 1. One of the goals 
of our work here is to test for consistency between the two experiments, by using a common 
framework to describe the SZ and foreground components. As this article was being prepared, 
refined spectra from SPT at 150 GHz were reported in [66]; we do not include these latest 
data in our comparison. 

Before fitting the SPT data with the ACT secondary model, we confirm that we recover 
the parameters and y 2 obtained for the model used in [54] using the SPT data. To combine 
the data over the full angular range, we follow the method in [54], using the [36] data at 
£ < 2000 (SPT-low) and the [54] data at smaller scales (SPT-high). More radio source power 
has been removed from the SPT-high spectra due to masking at a deeper flux level, so the 
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Figure 5. Distributions for secondary parameters from ACT and SPT, for best-fitting ACDM model. 
Parameters {a t sz, Oksz, a P , a c , a gs , a ge , a s j are the £(£ + l)Ce/2n power in /iK 2 at l = 3000 
and frequency 150 GHz. The tSZ-CIB correlation parameter £ is also defined at £ = 3000. The 
dust emissivity index /3 C is in flux units, for a modified blackbody with effective temperature 9.7 K. 
Conversions to power at each frequency are given in table 2. Strong priors, described in section 3.1.2, 
are imposed on {£, a ge , a gs , a s }. 


expected residual radio power in SPT-low is $3000 — 10.5 ± 2.4 /rK 2 , compared to 1.3 ± 0.2 
for SPT-liigh. We account for this by first subtracting a radio Poisson power of Be — 9.2 /iK 2 
from the SPT-low data, following the approach in [54], A Gaussian prior is then imposed on 
the overall residual radio level in SPT of a s / = 1.3 ± 0.2. 

We then extend the ACT secondary model to fit the SPT power spectra. Six of the ACT 
model parameters are expected to be common for the SPT data (the SZ and CIB parameters: 
atsz, flksz, £, Op, ac and /3 C ). In addition, to fit the SPT data we require a separate radio 
source parameter, a s ', and three calibration parameters, y\, j/ 2 , 2 / 3 , to calibrate the 95, 150, 
and 220 GHz maps respectively. We impose a uniform prior on these calibration parameters, 
as the SPT covariance matrices include the calibration uncertainty. 

The likelihood for ACT and SPT together is given by 


— 21n2z? = —2 ln^AC'r — 21n_£?sPT • 

The SPT likelihood is constructed as in eq. (3.1), with model and data vectors 


n — ^ 95,95 ^, 95,150 ^-, 95,220 ^, 150, 150 ^ 150, 220 ^ 220 , 220 , 

u b — A J 


(3.2) 

(3.3) 
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for SPT-liigh at l > 2000, and Cb = C^ 50,150 for SPT-low at t < 2000. To compute the 
model, we use the band-centers for SZ, radio and dusty sources given in table 1, from [54]. 

There is some degree of covariance between the ACT and SPT spectra, due to the 54 
deg 2 overlapping region of sky. Covariance between the two spectra due to cosmic variance, 
scaling as l// s ky> is estimated at a level of 8% (using 54/\/590 x 800); the addition of noise 
lowers this level, so we neglect the correlation in our combined analysis. 

3.3 Multi-frequency likelihood prescription 

To return the ACT (or ACT+SPT) multi- frequency likelihood for a given model we follow 
this approach: 

• Select primary cosmological parameters, and compute a theoretical lensed CMB power 

spectrum using the CAMB numerical Boltzmann code ([42]). 

• Select values for common secondary parameters: 0 — {a t sz> aksz> C a p , o, c , /3 C }. 

• Select values for ACT-specific secondary and calibration parameters: 9 — {a s , a ge , a gs , 
yi e , y 2 e, yis, U 2 s}, and/or SPT-specific secondary and calibration parameters: 9 — {a s /, 
Vi, V 2 , Vs}- 

• Compute the total theoretical secondary power spectra ^ ec,1J for all the required cross- 
spectra with eq. (2.4), using the effective frequencies for each experiment. 

• Compute the total model power at each frequency, Bf 1 '^ — Bf MB + ^ ec,lj . 

• Compute the bandpower theoretical power spectra for each dataset for both South and 
Equatorial regions for ACT (and for SPT), and compute the likelihood using eq. (3.1). 

3.3.1 Combining with WMAP 

In [61] we use the ACT and SPT likelihood in combination with data from WMAP seven- 
year data to estimate cosmological parameters. The WMAP data measure (. < 1000 angular 
scales, and so have minimal contamination from SZ and point sources. The public 7-year 
likelihood estimates the temperature spectrum from V and W bands (61 and 94 GHz, [40]). 
At these frequencies and angular scales, the infrared point source contribution is expected to 
be negligible, consistent with ACT and SPT measurements. The radio point source level is 
estimated and subtracted internally to the WMAP analysis using the multi-frequency data, 
e.g., [48]. Finally, in light of observations by both ACT and SPT, we also neglect the SZ 
power in the WMAP data, as it is expected to be small: 1(1 + l)Ciooo/27r ~ 12 yl\ 2 at 
61 GHz from the [4] model, assuming as = 0.8. 

4 Tests of the multi-frequency likelihood 

In this section we test the goodness of fit of the model to the ACT power spectra, assuming 
the ACDM cosmological model. We estimate the probability distributions of the secondary 
parameters using the MCMC method described in [16], fixing the ACDM parameters at 
best-fitting values. 6 * * We then investigate a set of possible extensions or modifications to the 

6 Physical baryon density = 0.02226, cold dark matter density f l c h 2 = 0.1122, ratio of the acoustic 

horizon to the angular diameter distance at decoupling 0 = 1.040, scalar amplitude ln[10 lo T s ] = 3.186, 

spectral index n s = 0.9707, and optical depth r = 0.898. 
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Figure 6. Power spectra measured by ACT at 148-218 GHz. with the best-fitting individual SZ and 
foreground components from table 2. The Galactic cirrus component has been subtracted. At 148 GHz 
(top) the secondary components are significant at scales smaller than l ~ 2000, with contributions 
from tSZ, kSZ, radio galaxies, the CIB, and the tSZ-CIB cross-correlation. The tSZ, kSZ and tSZ-CIB 
are non-zero in this model, but are not individually significantly detected from the ACT spectra. The 
radio power is constrained by bright source counts. At 218 GHz (bottom) the secondary signal is 
significant by l ~ 1000, and is dominated by the Poisson and clustered CIB. 


secondary model. We include the SPT power spectra and examine the consistency of the 
foreground model between the two datasets. We use the ACT ‘vl’ spectra for this analysis; 
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as described earlier, a refined estimate of the beams became available after the analysis was 
complete. We have checked that effects on parameters are negligible (< 0.2er), so do not 
update the parameter constraints or plots in this section. 


4.1 ACT data 

We find that the model provides a good fit to the ACT data over the full range of angular 
scales and frequencies. Figure 4 shows the total spectra (coadded over ACT-E and ACT- 
S, with the best-fitting Galactic cirrus component removed) decomposed into primary and 
secondary contributions. The SZ and foregrounds dominate at l ^ 2400 at 218 GHz, and 
at i ^ 3200 at 148 GHz. The goodness of fit is x 2 — 675 for 697 dof (reduced x 2 — 0.98, 
with PTE— 0.72, for 710 data points and 13 parameters). This indicates a good overall fit, 
but localised deviations can be hard to identify using the x 2 over the full angular range. 
Figure 4 therefore shows the residual power after subtracting the best-fitting model; we do 
not observe any significant features, indicating that the model fits both the angular and 
frequency dependence of the data in both regions. There is a positive excess in the ACT-E 
residuals at the smallest scales at 218 GHz, but this is consistent with correlated beam error, 
accounted for in the covariance matrix. 

The marginalized distributions for the secondary parameters fitting the data are shown 
in figure 5 and summarized in table 2. The Poisson-like and clustered CIB power, a p and 
a c , are detected at high significance, with index /3 = 2.2 ± 0.1 consistent with [2] who 
find 2.20 ± 0.07. The tSZ and kSZ power are individually seen at low significance, with 
an anti-correlation between atsz an d Uksz- The kSZ power peaks at a non-zero value, but 
the distribution is broad and consistent with zero. The total SZ power is detected at high 
significance. The tSZ-CIB correlation coefficient is unconstrained in the prior range 0 < £ < 
0.2, and is also unconstrained by ACT if allowed to vary over a broader range (e.g., £ < 0.5). 
The parameters for the power from radio sources and from Galactic cirrus are driven by their 
prior distributions. [61] present a physical interpretation of these parameters; the constraints 
are consistent with those found in the 1-year ACT analysis in [16], with reduced errors. 

In figure 6 we show the individual components that contribute to the 148 GHz and 
218 GHz power spectra after removal of the best-fitting Galactic cirrus power. At 148 GHz 
there are contributions from all the components. At 218 GHz the secondary spectrum is 
dominated by dusty point sources, both clustered and Poisson. This is illustrated further in 
figure 7, which shows the frequency dependence of the dominant components in our model 
at £ — 3000. The derived constraints on the CIB and radio source components, and the 
Galactic cirrus emission, at the ACT effective frequencies are also given in table 3 to allow 
comparison with other models. 


7 A flat prior is imposed, unless indicated as a Gaussian with x ± y for mean x and standard deviation y. 
s Results are reported as 68% confidence levels or 95% upper limits; £ is unconstrained so the prior upper 
limit is reported. 

9 If the prior on £ is broadened to 0 < £ < 0.5, the upper limit increases to dksz < 6.9 [61]. 

10 The SPT cirrus level we use is B 3000 = 0.16, 0.21, and 2.19 /rK 2 at 95, 150, and 220 GHz, as measured 
in [54]. 

11 Secondary parameters marginalized over the 6 ACDM model parameters are reported in table 1 of [61], 
and are consistent with these results. The marginalization has little effect on these secondary parameters, 
increasing errors by at most 10%. 

12 Gal-E and Gal-S are the Galactic cirrus powers in the ACT-E and ACT-S spectra. The levels are close 
to the priors imposed from the measured cross-correlations with IRIS ([14]). 
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Parameter 

Prior 7 

ACT 8 

SPT 

ACT±SPT 

sz 

atsz 

> 0 

3.3 ±1.4 

4.1 ±0.9 

4.0 ±0.9 


flkSZ 

> 0 

< 8.6 

< 4.2 

< 5.0 9 

CIB 

dp 

> 0 

6.9 ±0.4 

7.0 ±0.4 

7.0 ±0.3 


d c 

> 0 

4.9 ±0.9 

6.0 ±0.7 

5.7 ±0.6 


Pc 

> 0 

2.2 ±0.1 

2.0 ±0.1 

2.10 ±0.07 

tSZ-CIB 

£ 

0 < £ < 0.2 

< 0.2 

< 0.2 

< 0.2 

Radio 

d s 

2.9 ±0.4 

3.1 ±0.4 

— 

3.2 ±0.3 


d s ' 

1.3 ±0.2 

— 

1.4 ±0.1 

1.4 ±0.1 

Galactic cirrus 10 

a ge 

0.8 ±0.2 

0.9 ±0.2 

— 

0.9 ±0.2 


a gs 

0.4 ±0.2 

< 0.73 

— 

< 0.70 

Calibration 

Vie 

1.00 ±0.02 

1.010 ±0.007 

— 

1.006 ±0.006 


V2e 

— 

0.99 ±0.01 

— 

0.99 ±0.01 


Vis 

1.00 ±0.02 

1.011 ±0.007 

— 

1.010 ±0.007 


V2s 

— 

1.03 ±0.01 

— 

1.02 ±0.01 


yi 

— 

— 

1.01 ±0.02 

1.01 ±0.02 


V2 

— 

— 

1.007 ±0.008 

1.008 ±0.008 


V3 

— 

— 

1.02 ±0.02 

1.03 ±0.02 


best fit y 2 /dof 


675/697 

96/107 

773/810 


PTE 


0.72 

0.77 

0.82 


Table 2 . Likelihood parameters, assuming best-fit 6-parameter ACDM for the lensed CMB. 11 


ACT 

148 GHz 218 GHz 

95 GHz 

SPT 
150 GHz 

220GHz 

CIB-P 

CIB-C 

6.8 ±0.4 

4.8 ±0.9 

78 ±12 
54 ± 16 

0.90 ±0.02 
0.76 ±0.02 

8.0 ±0.5 
6.8 ±0.8 

69 ± 10 
59 ± 12 

Radio 

3.2 ±0.4 

1.4 ±0.2 

7.2 ±0.8 

1.4 ±0.2 

0.7 ±0.1 

Gal-E 12 

Gal-S 

0.9 ±0.2 
0.4 ±0.2 

11 ±2.3 
5.0 ±2.3 





Table 3. Derived constraints on foreground power, B3000 (a*K 2 ). 


4.2 Combination with SPT 

The same model also provides a good fit to the SPT spectra. The SPT data extend the 
frequency range to 95 GHz, adding three additional cross-spectra to the likelihood. We show 
the parameters estimated from SPT alone in figure 5 and in table 2; they are consistent with 
those from ACT, with ~ la shifts in aksz and j3 c . The radio Poisson level is lower due to the 
greater number of radio sources masked in the SPT maps. The model is shown with the SPT 
spectra in figure 8; the goodness of fit is \ 2 — 96 for 107 dof (reduced % 2 = 0.89, PTE=0.77). 
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Figure 7. Frequency dependence of the dominant components of the foreground power at l = 3000 
measured by the combined ACT and SPT data sets. The bands show the ler uncertainties from 
table 2. At 150 — 220 GHz the power from fluctuations in the CIB dominates; at lower frequencies the 
thermal SZ and radio source power is more significant. The SPT radio power is lower due to deeper 
integration. The kSZ and tSZ-CIB components are not shown. 


Given the consistency of the two datasets, we combine them to generate a joint likeli- 
hood; figure 5 includes the secondary parameters derived from a joint fit. In this case there 
are ten foreground parameters, and seven calibration parameters. The tenth foreground pa- 
rameter (not plotted) is a s ' for the Poisson radio sources in SPT. The goodness of fit of the 
joint model is \ 2 — 773, which can be compared to y 2 = 675 + 96 = 771 for the independent 
fit to each data set. This supports their consistency. 

We report the derived constraints on the CIB and radio source components at the ACT 
and SPT effective frequencies for each band in table 3. A difference of approximately 15% is 
expected between the CIB power at 148 GHz for ACT and 150 GHz for SPT, due to different 
effective bandpass frequencies and the strong CIB frequency dependence across the nun-wave 
bands. 

4.3 Tests of the likelihood 

This model fits the ACT and SPT data, and includes our uncertainties about the physical 
components, with priors describing our knowledge from other observations. However, it is 
a simplified parameterization of the emission. We therefore consider a set of extensions or 
modifications to the model, and test how the goodness of fit to the ACT data is affected 

13 Does not include calibration parameters. 

14 We report the Ay 2 to the nearest integer. 
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I i 

Figure 8. Power spectra at 95, 150 and 220 GHz, and their cross-spectra, measured by SPT ([54]), 
fit with the same model as the ACT data in figure 4. At 150 GHz the SPT-low spectrum from [36] is 
included, with excess radio power subtracted for comparison. Accounting for the different flux cuts 
applied to the ACT and SPT maps, and the different bandpass effective frequencies, the spectra are 
consistent. 


Model 

Number of 
parameters 14 

ACT 14 

A X 2 

Fiducial 

9 

0 

CIB index n free 

10 

-4 

Pc # Pp 

10 

0 

CIB Poisson corr m 0.8 

9 

5 

CIB Clustered corr = 0.8 

9 

1 

CIB T d = 13.6 K 

9 

1 

Fixed kSZ, aksz = 1-5 

8 

3 

Altered kSZ shape 

9 

1 

No tSZ-CIB corr, £ = 0 

8 

2 

No SZ 

6 

21 

Radio index a s — 0 

9 

1 

No Galactic residual 

7 

6 


Table 4. Modifications to secondary model. 


by an increase or decrease in parameters, or a change in the prior assumptions. In these 
tests, summarized in table 4, we hold the cosmological model fixed at the best-fitting ACDM 
parameters. A subset of these extensions are considered further in [61], testing their effect 
on the primary cosmological parameters. 

The CIB appears to be well-fit currently by a power-law in angular scale, with Bi oc 
£ 0 - 8 . [2] find an uncertainty of 0.06 in this scaling. If we allow the index n to vary, we find 
no improvement in the fit, but parameter distributions for the CIB parameters a p and a c 
are broadened, as they are correlated with the power-law scaling. We also assume that the 
CIB emission is perfectly correlated among frequencies, in the range 95-220 GHz. Evidence 
for imperfect correlation was reported in [52]. The effect of this assumption is tested by 
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setting the correlation coefficient to < 1 in the model, choosing 0.8 for either the Poisson 
or clustered components, roughly corresponding to the degree of correlation between maps 
reported in [52]. For the Poisson component we find that this degrades the goodness of fit 
by Ay 2 = 5 compared to the perfectly correlated case. We also assume a common frequency 
scaling of the clustered and Poisson terms. A different scaling in frequency may be expected 
if, for instance, the redshift dependence of the clustered and Poisson power is different, so 
a common index is not necessarily expected. Allowing it to vary independently does not 
significantly improve the goodness of fit, but does lead to a poorly constrained distribution 
for /3 C , and removes the detection of the CIB at 148 GHz. The index for the Poisson sources 
in this case, /3 P — 2.14 ± 0.15, is consistent with the joint index. Changing the effective dust 
temperature from 9.7 K to 14 K, consistent with the value obtained by [24] from a fit to 
the FIRAS CIB frequency spectrum, has no effect on the model, apart from a corresponding 
change in /3 C . 

Exploring the SZ assumptions, we consider fixing the kSZ contribution to the linear 
theory estimate for a universe with as — 0.8, assuming homogeneous reionization. This 
degrades the goodness of fit of the model by only Ay 2 = 3, indicating that the data cannot 
yet distinguish between homogenous and patchy reionization. Limiting the kSZ in this way 
also leads to tighter constraints on the tSZ power. Modifying instead the shape, we find that 
adding a patchy reionization template from [6] to the kSZ, which changes its shape, does not 
affect the other secondary parameters. Neglecting the tSZ-CIB cross-correlation also fits the 
data equally well for one fewer parameter. The dependence of the kSZ constraints on the 
model for the tSZ-CIB correlation is explored further in [61]. If we neglect the SZ components 
altogether, setting atsz m «kSZ — 0 and keeping only the CIB and Galactic components, the 
goodness of fit significantly worsens, with an increase of Ay 2 = 21. 

Our model imposes an a priori assumption on the frequency scaling of the extragalactic 
radio sources. We test the effect of changing the radio spectral index to a s — 0, finding 
negligible effect on parameters and goodness of fit. There is also little effect from changing 
the prior on the power from [26] by lcr, corresponding to a different model for the bright 
radio sources that lie above the detection threshold. Removing the prior altogether opens up 
degeneracies with other parameters, but does not significantly improve the goodness of fit. 

Finally, we test the effect of removing the Galactic cirrus components; the goodness of 
fit worsens by Ay 2 = 6 and the clustered CIB level increases. This indicates a preference for 
Galactic cirrus at the 95% confidence level. 

5 CMB-only likelihood 

Understanding the contribution of the secondary components to the ACT power spectra is 
vital for extracting the cosmological information, due to possible degeneracies between pri- 
mary and secondary parameters. Values of the secondary parameters are also astrophysically 
interesting. However, if we are only interested in the cosmological parameters, a simplified 
likelihood is desirable. 

We construct a CMB-only likelihood from ACT data as follows. Instead of using the 
ACT likelihood to estimate cosmological parameters, we take the intermediate step of esti- 
mating the CMB power spectrum in bandpowers, marginalizing over the possible contam- 
ination. This is a natural extension to forms of CMB data compression that have been 
adopted in earlier analyses, e.g., [10]. Such ‘grand unified spectra’ were used in a number 
of subsequent papers to combine the results from various CMB experiments, marginalizing 
over a variety of nuisance parameters, e.g., in [8, 9, 62], At large scales, where contamination 
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Figure 9. Estimated CMB bandpowers from ACT, marginalized over extragalactic source and SZ 
components. Bandpowers are estimated for ACT-E and ACT-S separately; here we show the inverse- 
variance weighted combination. The bandpowers are correlated at the ~ 20% level at scales l £ 2000 
due to covariance with the secondary parameters. The total multi-frequency spectra for ACT-E 
(dashed, at 148 GHz, 148x218 GHz, and 218 GHz) are also shown to indicate the significant level of 
SZ and foreground power at small scales. 

from SZ and point sources is negligible, the estimated CMB is simply an optimally combined 
average of the multi- frequency spectra as in e.g., [32]. At smaller scales the CMB spectrum 
has additional uncertainty due to secondary contamination. 

By marginalizing over nuisance parameters in the spectrum-estimation step, we can 
effectively decouple the primary CMB from non-CMB information. No additional nuisance 
parameters are then needed when estimating cosmological parameters. 

5.1 Method: bandpowers via Gibbs sampling 

To implement this method in practice, we estimate n\, CMB bandpowers, marginalizing over 
the secondary parameters. We use the full multi-frequency likelihood from section 3 to do 
this, but estimate CMB bandpowers instead of cosmological parameters. 

We recall that the model for the theoretical power for a single cross- frequency, cross- 
season spectrum, , is written as 



where C s ^ c ' l \9) is the secondary signal as in eq. (2.4), and is a function of secondary param- 
eters 9. Writing the spectrum in bandpowers, where w are the bandpower 

window functions, we write the model for the bandpowers in vector form as 
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where C b h and C b ec are multi- frequency, multi-season vectors of length n b x u spec , where rib 
is the number of bandpowers, and n spe c is the number of cross-season and cross-frequency 
spectra (n spe c — 9 for ACT-E, and 18 for ACT-S). The secondary spectra differ between 
frequencies but not between seasons. The mapping matrix A, with elements that are either 
1 or 0, maps the CMB bandpower vector (of length rib), which is the same at all frequencies 
and in all seasons, onto the ( rib x ^s P ec)-length data vector. 

We want to estimate CF MB , marginalized over the secondary parameters, 9. The poste- 
rior distribution for C)p MB , given the observed multi- frequency, multi-season spectra Cb, can 
be written as 

p(C£ MB \C b ) = J p(C£ MB ,9\C b )p(9)d9. (5.3) 

We find that Gibbs sampling provides an efficient way to map out the joint distribution 
p(C b MB , 9\Cb), and to extract the desired marginalized distribution p(C'^ MB |C; ) ). 

Gibbs sampling can be used in the special case that at least one conditional slice through 
a multi-dimensional distribution has a known form, and has been used, for example, to 
estimate the large-scale CMB power spectrum, and to marginalize over Galactic foregrounds, 
e.g., [17, 20, 35, 40, 72]. Here, we split the joint distribution into two conditional distributions: 
p(C'^ MB |0, Cb), and p(9\C b !MB , Cb). We write the multi-frequency likelihood for a single ACT 
region, from eq. (3.1), as 

- 21nJ^ = (AC™ b + Cf c - C fe ) T S- 1 (AC fe GMB + C s b ec - C b ) 

+ lndetS, (5.4) 

which is a multivariate Gaussian. 

If C b ec is held fixed, the conditional distribution for the CMB bandpowers, 
p(C b MB \9,C b ), assuming a uniform prior for p(C)p MB ), is then also a Gaussian. It has a 
distribution given by 

- 21np(C£ MB \9,C b ) = (C fe CMB - C 6 ) T Q" 1 (C' b CMB - C b ) 

-blndetQ, (5.5) 

The mean, G b , and covariance, Q, of this conditional distribution are obtained by taking the 
derivatives of the likelihood in eq. (5.4) with respect to cF MB . This gives mean 

Cb = [A t E- 1 A]- 1 [A t E- 1 (C 6 - Cf c )], (5.6) 

and covariance 

Q = A t S“ 1 A. (5.7) 

We can draw a random sample from this Gaussian distribution by taking the Cholesky 
decomposition of the covariance matrix, Q m LL r , and drawing a vector of Gaussian random 
variates G. The sample is then given by (3)p MB = C b + L G . 

If instead C)p MB is held fixed, the conditional distribution for the secondary parameters, 
p(9\C™ B ,C b ), is not a Gaussian, but can be sampled with the Metropolis algorithm that 
is used in the MCMC sampling in section 4. To map out the full joint distribution for 9 
and C fc GMB we alternate a Gibbs sampling step, drawing a new vector of CMB bandpowers, 
with a Metropolis step, drawing a trial vector of the secondary parameters 9. 

We choose a uniform positive prior distribution for p(CF MB ), and restrict the CMB 
bandpowers to be zero at £ > 4500, where the CMB power is expected to be less than 1 pK 2 . 
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About 100,000 steps are required for convergence of the joint distribution, assessed with 
the [15] spectral test. The mean and covariance of the resulting marginalized bandpowers, 
C^ MB , are then estimated following the standard MCMC prescription, e.g., [41, 65]. 

5.1.1 Combining spectra from different regions 

There is only one underlying CMB power spectrum, so this method could be used to estimate 
a single spectrum, or set of bandpowers, from the two ACT regions. However, the bandpower 
window functions are different for each region due to their distinct geometries. To easily con- 
serve this information, we estimate the CMB bandpowers for ACT-E and ACT-S separately. 
Since the secondary parameters are common to both, the estimated CMB bandpowers will 
be correlated between the regions at small scales. 

To estimate the joint distribution for the ACT-E and ACT-S bandpowers, we map out 
p(C™ B - E ,C™ B ~ S ,e\C b ) by taking sequential sampling steps from the conditional distri- 
butions: 

p(c fe CMB - E |c b CMB - s ,0,c fe ), 

p{c CMB-S\ c CMB-E^c b)j 

p(e\C™ B - E ,C™ B ~ s ,C b ). (5.8) 

The marginalized distribution for the CMB bandpowers, p(C^ MB ~ E , C^ MB-S |Cfc), with its 
associated covariance matrix, is then computed from the samples. This could be extended 
to include the SPT data, or data from Planck, for example. 

5.1.2 Calibration factors 

There are four ACT calibration factors. To minimize bin-to-bin correlations in the estimated 
CMB bandpowers due to calibration uncertainty, we divide out the 148 GHz calibrations for 

the two ACT spectra, estimating C( CMB ~ E = C^ MB ~ E /yf e for the ACT-E bandpowers, and 
c /CMB-S = qCMB— S / y 2 s for ACT _ S . 

We then estimate the 148 GHz calibration factors, y\ e , y\ s , and relative 218/148 GHz 
calibration factors, y 2 e/yie, V 2 s/yisi as part of the secondary parameter set. 

5.2 Marginalized CMB bandpowers 

Figure 9 shows the estimated CMB bandpowers from the ACT-E and ACT-S spectra, co- 
added together and compared to the multi-frequency spectra. The bandpowers for each 
region are reported in table 5. In this table we report the CMB spectra derived using the 
updated ‘v2’ multi-frequency spectra. Without assuming any cosmological model, the CMB 
bandpowers over the full angular range are remarkably consistent with the theoretical ACDM 
model predicted by WMAP. The uncertainty on the bandpowers rises at scales smaller than 
£ ~ 3000, and the correlations between bandpowers increases. 

Figure 10 shows the effect of marginalization on the bandpower errors, using the ratio 
between the marginalized errors and the unmarginalized errors for a fixed secondary model. It 
is clear that by measuring the spectrum at multiple frequencies, the CMB can be successfully 
separated from secondary contamination out to scales £ ;$ 3500. At scales £ < 2000 there is 
little error inflation due to foreground uncertainty, and the errors are inflated by ^20% (15%) 
by £ — 3200 for the ACT-E (ACT-S) spectra. The marginalized distributions for the CMB 
bandpowers are well approximated by Gaussians for multipoles to band-center £ — 3540 for 
ACT, as shown in the appendix. 


22 


J CAP 0 7 (2013) 025 


4 

l{l+l)C h /2ir (/rK 2 ) 

ACT-E ACT-S Coadd 15 

590 

2157 ± 159 

2343 ± 160 

2250 ± 113 

690 

1729 ± 115 

1744 ± 107 

1737 ± 78 

790 

2499 ± 146 

2274 ± 126 

2370 ± 96 

890 

1945 ± 109 

1903 ± 102 

1923 ± 74 

990 

1068 ± 58 

1187 ±61 

1124 ±42 

1090 

1206 ± 61 

1149 ±58 

1176 ±42 

1190 

1036 ±51 

1016 ±48 

1026 ± 35 

1290 

679 ± 33 

766 ± 38 

717 ±25 

1390 

819 ±39 

787 ± 36 

802 ± 26 

1490 

661 ± 31 

650 ± 29 

655 ± 21 

1590 

452 ± 19 

474 ± 21 

462 ± 14 

1690 

387 ± 16 

355 ± 17 

372 ± 11 

1790 

344 ± 14 

347 ± 16 

345 ± 10 

1890 

242 ± 10 

246 ± 12 

244 ±8 

1990 

230 ± 10 

214 ± 11 

223 ±8 

2090 

199 ±9 

204 ± 11 

201 ± 7 

2240 

137 ±5 

127 ±6 

133 ±4 

2440 

92.9 ±3.9 

87.8 ±5.1 

91.0 ± 3.1 

2640 

57.6 ±3.3 

56.4 ±4.5 

57.2 ± 2.6 

2840 

43.0 ±3.2 

44.7 ±4.2 

43.6 ± 2.5 

3140 

22.5 ±2.2 

19.3 ±2.8 

21.3 ± 1.7 

3540 

9.2 ±2.3 

9.0 ±2.8 

9.1 ± 1.8 


Table 5. Lensed CMB anisotropy power. 10 


We compare the secondary parameters recovered in this model-independent sampling to 
the case where ACDM is assumed. This comparison is shown in the appendix; the parameters 
are consistent, with about a la shift in the estimated kSZ power. We find that the CMB 
bandpowers are not strongly correlated with the secondary parameters until scales well into 
the Silk damping tail at £ xi 2500; a dominant correlation is then with the kSZ power due to its 
blackbody frequency dependence, and a smaller kSZ power — compensated by larger primary 
CMB power — is allowed when the ACDM assumption is relaxed. The CMB bandpower 
covariance matrix conserves this correlation information. 

5.3 The CMB-only likelihood 

We construct the CMB-only likelihood from the angular range where the CMB bandpowers 
are Gaussian, conservatively choosing i < 3500. We do not use the 3500 < i < 4500 
bandpowers as they are increasingly non-Gaussian, due to the foreground marginalization, 
and are more strongly correlated with foreground parameters. The likelihood is given by 

-21n^f((7 6 CMB |Q th ) = a: T 5r 1 x + lndetS. (5.9) 

ls This coadds the ACT-E and ACT-S CMB bandpowers for plotting purposes. 

16 To compute a likelihood using these data, ACT-E and ACT-S should be used with the covariance matrix 
and bandpower window functions provided on LAMBDA. 
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Figure 10. Inflation of errors due to foreground marginalization, relative to the errors for a best- 
fitting foreground model: at scales smaller than l ~ 2000, the errors are increased due to foreground 
uncertainty. 


Here 


_ ( C™ B ~ E - U> M ,ACT-E Cf \ 
V C fe GMB - S - WU.ACT- S cf ) 


(5.10) 


where C]p MB anc [ ^ are the marginalized mean and covariance matrix for the bandpowers, 
and Cf 1 is the lensed CMB spectrum generated from e.g., CAMB. We use 21 bandpowers 
for ACT-E and ACT-S in the range 500 < £ < 3500. A single calibration parameter for 
each region is marginalized over analytically, following [11], The prescription for using this 
likelihood is simple, as no extra nuisance parameters are needed. 

To test the performance of this compressed likelihood, results are compared using the 
full multi-frequency likelihood, and the CMB-only likelihood. Cosmological parameters are 
estimated for the restricted ACDM 6-parameter model, and a set of more extended models 
that probe the damping tail and peak shapes, including the running of the spectral index, the 
number of relativistic degrees of freedom N e g, the lensing amplitude Al, and the variation in 
fine structure constant, a. Parameter constraints using both likelihoods agree to O.lcr, and 
are reported in [61]. We conclude that this is an efficient alternative to the full likelihood for 
the typical extensions considered in cosmological analyses, although the full likelihood may 
give more optimal results for unusual models with features far into the damping tail. 


6 Summary 

In this paper we have presented a likelihood formalism to describe the ACT multi-frequency 
power spectra that includes contributions from SZ and foreground components in addition 
to the lensed CMB. We model the data including four late-time astrophysical components: 
thermal and kinetic SZ, emission from CIB galaxies, and emission from radio galaxies. 

We have quantified these components using seven power spectra, splitting the CIB 
into a Poisson and clustered part, and including power from the cross-correlation between 
tSZ emission from clusters, and emission from CIB galaxies that also trace the large scale 
structure. Rather than a minimal model with the fewest parameters, we have sought a 
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model that includes our uncertainties with priors describing our knowledge from additional 
observations. For example, while the data do not demand that we include the tSZ-CIB 
correlation, we are motivated to include it to avoid placing unphysically strong limits on the 
kSZ power. 

Modeling these astrophysical components allows us to probe the primordial CMB fluc- 
tuations down to an angular resolution of 4’ using ACT. We have used the model to extract 
an estimate of the primordial CMB spectrum well into the Silk damping tail, marginalizing 
over the foreground uncertainty. This produces a simplified compressed likelihood for use in 
cosmological parameter estimation. 

We find that data observed by the South Pole Telescope give results consistent with 
ACT, accounting for the different removal of radio point sources, and different degree of con- 
tamination by Galactic cirrus. SPT and ACT have very different instrument design and scan 
strategies, and their observations on the sky have limited overlap. The excellent agreement 
between the datasets is not only an important cross-check but is another demonstration of 
cosmic homogeneity. 
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A Tests of the CMB likelihood 

In this appendix, we perform additional tests on the CMB-only likelihood. In figure 11 we 
show a selection of the distributions of the CMB bandpowers from the estimated 600 < £ < 
4500 range. Distributions for the ACT-E and ACT-S bandpowers are compared to Gaussian 
distributions (dashed curves); bandpowers at £ > 3900 are significantly non-Gaussian, but 
are well fit by Gaussians at larger scales. The same behaviour is found for the ACT-S 
bandpowers. 

We then compare the secondary parameters estimated in two ways: (1) estimating 
CMB bandpowers, and (2) estimating 6 ACDM parameters. The distributions are shown 
in figure 12, and are consistent. The tSZ, point source parameters, and Galactic cirrus 
parameters are not strongly affected by the CMB model assumptions. The kSZ power, aksz, is 
la lower in the model-independent case, as it is anti-correlated with the CMB bandpowers 
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Figure 11. Probability distributions of a selection of CMB bandpowers for ACT-E in the range 
600 < £ < 4500, marginalized over secondary parameters (solid). The bin-center for each bandpower 
shown is indicated on each panel. The bandpowers are well-approximated by Gaussian distributions 
(dashed), except at scales £ 3700. The same behavior is seen for the ACT-S bandpowers. 





Figure 12. Secondary SZ and foreground parameters estimated assuming the ACDM model 
(‘LCDM’), compared to the same parameters estimated jointly with primary CMB bandpowers (‘C’). 
They are consistent, but an anti-correlation between the primary CMB bandpowers and the kSZ 
power leads to a reduction in a^gz in the latter case. 


at £ > 2000 due to the common blackbody dependence. The data cannot distinguish between 
lensed CMB power and kSZ power at £ ~ 3000 scales, so the preference for a smaller kSZ 
value in the model-indepedent case is driven by the prior that the CMB power is positive. 
We test this by allowing the CMB bandpowers to take unphysical negative values. Here, the 
kSZ power increases to aksz < 12 at 95% confidence, more consistent with the limits when 
ACDM is assumed. 
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